Міністерство освіти та науки України Національний університет “Львівська політехніка” кафедра ЕМА Розрахункова робота “Математичне моделювання електромеханічних пристроїв”  Математична модель подвоювача частоти 1. Завдання. Скласти математичну модель електромагнітних перехідних процесів для подвоювала частоти, схема якого зображена на рис. 1
 рис. 1. Подавач частоти Характеристика намагнічування сталі задана: Вебер-амперна характеристика дроселя задана: Напруги живлення є задані функціями часу: Вхідні дані: Сумарні опори обмоток підмагнічування та первинних обмоток: Сумарні індуктивності розсіяння обмоток підмагнічування та первинних обмоток: Кількість витків обмоток підмагнічування, первинних та вторинних обмоток: Довжина силової лінії елементів магнітопровода: Площа поперечного перетину магнітопроводів: Розрахункова робота з предмету "Математичне моделювання електромеханічних перетворювачів | Виконав: студент групи ЕМА-41 Осадця О.П. !-------------------- Головний модуль програми ------------------------| ! Підключення модуля даних Use EMA ! Описова частина Dimension X_POPR(5),Y(4) ! Читання вхідних даних, початкової умови та вивід їх у файл вихідної ! інформації CALL DANI ! Час початку інтегрування T=T0 ! Запис у файл вихідної інформації даних сформованих на основі ! початкової умови CALL REZ(T) ! Формування кореня 0-го кроку інтегрування CALL POPER(X_POPR) ! Початкове значення лічильника кроку введення результатів LIC=0 ! Мітка початку кроку інтегрування 100 CONTINUE ! Визначення часу на кроці інтегрування T=T+h Write(*,*)' Час на кроці інтегрування T= ',T ! Перевірка умови: якщо час на кроці інтегрування більший за час ! припинення інтегрування, то слід припинити роботу програми : IF(T.GT.Tp) GOTO 110 ! Визначення значення напруг на кроці інтегрування CALL VOLT(T) ! Занулення лічильника ітерацій NI=0 ! Мітка початку ітерації методу Ньютона 200 CONTINUE ! Визначення номеру ітерації NI=NI+1 ! Перевірка умови: якщо кількість ітерацій перевищує 25 (це означає, ! що ітераційна процедура розбігається), слід зробити аварійну ! зупинку програми IF(NI.GT.25) GOTO 120 ! Обчислення нев'язки CALL NVZ(X_POPR,Y) ! Обчислення похідної нев'язки CALL MATR_A ! Визначення поправки первинних невідомих CALL GAUSS(A,Y,4,4,DELTA_XP) ! Визначення первинних невідомих CALL PR_NV ! Визначення вторинних невідомих CALL FORMUL ! Прийняття рішення про припинення ітераційного процесу call NSTOP(XP,DELTA_XP,4,NI,KZ) ! Перевірка умови: якщо KZ=0, то виконуємо наступну ітерацію IF (KZ.EQ.0) GOTO 200 write(*,*)'ni',ni WRITE(*,*)'Кількість ітерацій на кроці інтегрування NI= ',NI ! Збільшення значення лічильника кроку виведення результатів LIC=LIC+1 ! Перевірка умови: якщо значення лічильника дорівнює кроку виведення ! результатів, то вони записуються у файл IF (LIC.EQ.krok) THEN ! Запис результатів обчислень CALL REZ(T) ! Занулення значення лічильника кроку виведення результатів LIC=0 END IF ! Формування кореня кроку інтегрування CALL POPER(X_POPR) ! Повернення на початок наступного кроку інтегрування GOTO 100 ! Мітка аварійної зупинки програми 120 Write(*,*)' Аварійна зупинка програми NI>25 ' ! Закриваємо файл результатів CLOSE(20) STOP ! Мітка припинення інтегрування 110 Write(*,*)' Розрахунок закінчено ! ' ! Закриваємо файл результатів CLOSE(20) END   !_____________________________________________________________________| ! Підпрограма для читання вхідних даних, початкової умови та | ! виводу їх у файл вихідної інформації | !_____________________________________________________________________| SUBROUTINE DANI !_____________________________________________________________________| ! Вхідні дані: немає | ! Вихідні дані: | ! XP0 - вектор початкових значень первинних невідомих | ! XV0 - вектор початкових значень вторинних невідомих | !_____________________________________________________________________| ! Підключення модул...